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ancilla 
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Electron and nuclear spins associated with point defects in insulators are promising systems for solid 
state quantum technology 1 " 3 . While the electron spin usually is used for readout and addressing, 
nuclear spins are exquisite quantum bits 4 5 and memory systems 3 ' 6 . With these systems single-shot 
readout of nearby nuclear spins 5 as well as entanglement 4,7,8 aided by the electron spin has been 
shown. While the electron spin in this example is essential for readout it usually limits nuclear spin 
coherence 9 . This has set of the quest for defects with spin-free ground states 810 . Here, we isolate a 
hitherto unidentified defect in diamond and use it at room temperature to demonstrate optical spin 
polarization and readout with exceptionally high contrast (up to 45%), coherent manipulation of an 
individual excited triplet state spin, and coherent nuclear spin manipulation using the triplet electron 
spin as a meta-stable ancilla. By this we demonstrate nuclear magnetic resonance and Rabi 
oscillations of the uncoupled nuclear spin in the spin-free electronic ground state. Our study 
demonstrates that nuclei coupled to single metastable electron spins are useful quantum systems 
with long memory times despite electronic relaxation processes. 

Coupled electron and nuclear spins in solids are promising candidates for quantum memories and 
registers and present a particular class in the emerging field of hybrid quantum systems 711-14 , in which 
different types of qubits perform distinct functions according to their advantageous properties. In this 
particular case, the nuclear spin, which is weakly coupled to the environment, serves as a long-lived 



quantum memory, whereas the electron spin, which has a short coherence time, but interacts strongly 
with external fields, serves as a readout gate 4,7,9 . In this architecture, the permanent presence of the 
electron spin is a source of decoherence to the nuclear spin 8 . To avoid such electron spin induced 
decoherence, the electron that carries the spin should be physically removed when it is not needed and 
returned to the nuclear spin only when the initialization or readout gate is applied. One approach to 
address this problem is to frequently remove the electron, e.g. by photo-ionization with a strong laser, 
resulting in a motional narrowing type of decoupling 9 . Another approach - that is followed in this work - 
is to choose an electronic system with a spin-free ground state and a meta-stable excited spin state that 
can be optically pumped and activated prior to the application of the readout gate 8 . Such a system 
implemented with a single solid state defect would enable a universal nuclear spin gate. In the following 
we present a suitable defect in diamond and demonstrate coherent control of a single nuclear spin in 
the spin-free electronic ground state. Notably, besides the NV center, this is the second known single 
spin defect in the solid state, whose coherent spin motion is detectable at room temperature. 

The defect 1 is observed in ultra-pure single crystalline HPHT diamond after fabrication of vertical nano- 
wires (see Fig. la,b), whose main purpose is to provide high photon collection efficiency 15 . The wires are 
etched into the sample surface by reactive ion etching (see SI). Fig. lc shows a confocal fluorescence 
microscope image of the sample. To identify candidates for meta-stable triplet spin defects, we observe 
the fluorescence intensity of bright spots on the sample surface with and without a strong applied 
magnetic field from a permanent magnet (=0.5 T) (see Fig. Id). An arbitrary oriented magnetic field 
causes spin mixing and thereby modifies the steady-state fluorescence 16 . We observe two types of 
emitters with decreasing and increasing fluorescence response, corresponding to NV" centers and the 
new ST1 defects, respectively. The increase of the steady-state fluorescence is due to the mixing of 
meta-stable triplet states as will be shown later. 

To characterize the defect further, we record the photon auto-correlation and observe a dip well below 
0.5 at zero delay, which confirms that we indeed observe single defects 17 (see Fig. lg). At higher 
excitation intensity, the photon auto-correlation shows pronounced 'bunching-shoulders', which 
supports the idea of the existence of a long-lived shelving state. To characterize the defect fluorescence, 
we present fluorescence emission spectra at room-temperature and at 10K (see Fig. le). In both cases a 
strong resonant emission line is observed at 550(1) nm, with a broad phonon side band extending to the 
near infrared. Notably, this resonant transition does not correspond to a known diamond defect line 18 . 
To characterize the strength of the optical transition, we present the excited state lifetime that we have 
determined in a room temperature experiment by pulsed optical excitation (see Fig. If). An excited state 
lifetime of 9.5(2) ns is observed. 

Next, we search for a signature of an electron spin resonance in the optical fluorescence. Fig. 2a shows 
the fluorescence intensity that is observed while sweeping the microwaves from 200 to 1400 MHz. 
Three distinct lines marking an increase in fluorescence intensity upon spin resonance are observed at 
frequencies 2E = 278(1) MHz, D-E=996(l) MHz, D+E=1274(l) MHz. This optically detected magnetic 
resonance (ODMR) spectrum indicates a triplet electron spin with zero field splitting parameters 19 



1 hereafter termed STl-defect following the convention that is commonly applied in the EPR community 



D=1135(l) MHz and E=139(l) MHz. The observed large E parameter could originate from the presence 
of a large local strain or from an intrinsic low symmetry of the defect. We have compared the spin 
resonance spectra of ~20 individual defects (see SI) and observed variations of the E parameter smaller 
than 5 MHz, indicating that the E originates from an intrinsic low symmetry of the defect. From ODMR 
spectra with a variable axial magnetic field applied (Fig. 2b), we determine an electronic g-factor 
g e = 2.0(1), close to the ^-factor of the free electron spin. The absence of an intrinsic hyperfine 
splitting indicates that the defect's constituent atoms possess no nuclear spin (see SI for a detailed 
discussion of the possible constituents and structure of the defect). The aligned axial magnetic field has 
an angle of 35(1)° w.r.t the [111] crystal direction, indicating that the defect's main symmetry axis is 
oriented along the [110] crystal direction (see inset in Fig. 2b). 

All findings reported up to here are compatible with two possible electronic configurations: the triplet 
state could either be the long lived shelving state or the ground state. To show that the triplet state is 
indeed the long lived shelving state, we measure the decay of the shelving state population. If the 
shelving state was a singlet, a single exponential decay would be expected. By contrast, in the present 
case, three distinct exponential decay terms are observed (Fig. 3b) and their amplitudes are sensitive to 
microwave fields resonant with the triplet ESR transitions. In each decay curve shown, the dominant 
decay terms are controlled by switching the spin population with microwave n-pulses resonant with the 
three ESR transitions. This indicates that the shelving state is indeed the triplet state. The result can be 
used to assign the decay constants to the spin quantum numbers (Fig. 3c). A complete model of the 
optical cycle can be established by determining also the triplet population rate by an analogous pulsed 
experiment (see SI). Fig. 3d summarizes all results related to the triplet kinetics. The essential features 
are: a rapid (~l/70 ns) inter-system-crossing rate from the singlet excited state to the triplet state and 
substantial differences among the triplet lifetimes, Tio^SOO ns, T|_i>=1000 ns, T|+i>=200 ns. These are 
responsible for the positive and exceptionally high (45 %) optical spin signal contrast, and increasing 
steady-state fluorescence response to a randomly aligned static magnetic field as observed in Fig. Id 
because mixing the triplet sublevels results in faster decay to the ground state. It follows that optical 
pumping prepares the system in a spin polarized state of the meta-stable triplet state with predominant 
population of the slowly decaying triplet sublevel, 1 0). We now apply coherent spin manipulation to 
demonstrate Rabi oscillations between the triplet spin states (Fig. 4). Because the spin population of the 
faster decaying states | +) or | -) is increased by a resonant ESR pulse, any coherent spin manipulation 
results in the enhancement of the fluorescence intensity as shown in Fig.4b. The oscillations are 
described by two exponential decay functions, one of which is modulated by a cosine. The decay times 
and oscillation amplitudes correspond to the average of the lifetimes and the difference of the initial 
populations of the two involved spin states, respectively, and are in excellent agreement with the rates 
measured independently above (Fig. 3d). Note that the decay times are fully accounted for by the triplet 
lifetimes, which implies that spin relaxation is negligible. The triplet lifetimes are long enough to 
coherently address a coupled nuclear spin. We have compared the Zeeman sweeps of about 20 defects 
and observe that two of them show hyperfine splittings due to proximal / = 1/2 carbon nuclear spins of 
several tens of MHz see (Fig. 5a). 



We now show the application of the metastable electron spin as a universal single nuclear spin gate and 
demonstrate coherent manipulation and readout of the nuclear spin in the spin-free electronic ground 
state. To achieve this, we exploit the metastable electronic triplet spin ancilla to transfer spin- 
polarization onto and off the nuclear spin. To this end we exploit the optical nuclear polarization (ONP) 
induced by a level-anti-crossing (LAC) between the | -l,+l/2) and 1 0,-1/2) hyperfine levels at about 42 
mT to selectively mix electronic and nuclear spin states 2020 (see SI). Under these conditions, spin- 
polarization is transferred to the nuclear spin when the electron enters the meta-stable triplet level 
resulting in optical polarization of the nuclear spin 21 . After the optical pump pulse is turned off, the 
electron decays to the singlet ground state within a time given by the triplet lifetime. This decay 
eliminates the triplet electron spin and leaves the nuclear spin in its polarized state. Likewise, when the 
system is optically pumped from the singlet ground state to the excited triplet state, the nuclear spin 
state is preserved and the hyperfine levels are populated correspondingly. At the LAC, the different 
lifetimes of the hyperfine levels result in different fluorescent intensity of the nuclear spin projections. 
This provides a mechanism to readout the nuclear spin (see Fig. 5b and Methods). Fig. 5c shows the 
nuclear spin signal contrast as the energy levels are tuned into the LAC. Nuclear spin polarization is 
maximal at the LAC and vanishes when the detuning from the LAC becomes comparable to the hyperfine 
coupling strength. As an application of the nuclear spin manipulation, we show nuclear magnetic 
resonance (NMR) and coherent manipulation of the nuclear spin via Rabi oscillations (Fig. 5 d,e). For the 
measurement of a nuclear magnetic resonance spectrum, we apply a radio-frequency (RF) n-pulse of 
variable frequency. We observe a resonance dip at 450(1) kHz, corresponding to the Larmor precession 
frequency v 13C = 7i 3C 5 = 447 5 kHz of a bare 13 C nucleus. Here, y 13C = 10.705 kHz/mT is the 
gyromagnetic ratio of the 13 C nuclear spin and B = 41.8(5) mT is the applied field (corresponding to the 
LAC). We can thereby identify the coupled nuclear spin as a 13 C lattice spin. Note the absence of a 
chemical shift which indicates that indeed the electron spin is absent. For the measurement of nuclear 
Rabi oscillations, we apply a resonant RF pulse of variable length between the initialization and readout 
laser pulses. Several Rabi oscillations are observed that decay on a time scale of about 1 ms. The decay 
is likely to originate from dipole-dipole coupling of the nuclear spin to surrounding 13 C lattice spins. Note 
that the dominant decay mechanism is spin diffusion. By contrast to ensemble magnetization (that 
decays with TJ, here the relevant time scale is the ensemble T 2 * time that is about 10 ms in diamond 
with natural 13 C abundance 22 and agrees reasonably with the observed decoherence time. For future 
increase of the coherence time, such dipole interactions between nuclei of the same species can be 
efficiently removed by homonuclear decoupling schemes 9,23-25 . Other sources of decoherence could be 
paramagnetic or nuclear spin impurities at the surface whose decoherence effects could be avoided by 
isolating ST1 defects in the bulk crystal. 

The defect spin system presented in this work holds promise for room temperature solid state nuclear 
spin memories and quantum registers with ultra-long spin coherence times. Such a device could be used 
as a memory in integrated hybrid quantum systems 1314 or as a solid state quantum token 9 . For such 
applications, coherence times on the order of seconds or hours are desirable. Owing to its spin-free 
ground state, the system presented here is a promising candidate for achieving such long coherence 
times by homonuclear dynamical decoupling, such as BLEW 23 , FSLG 24 , MREV 25 . An increase of the 
nuclear spin coherence time by several orders of magnitude is expected 9 . In diamond, the coherence 



time of 13 C nuclear spins (at room temperature) is ultimately limited to time scales on the order of one 
day 22 . 



Methods 

Polarization and readout of the triplet spin state 

A long intense laser pulse is applied to prepare the system in a spin-polarized state of the triplet level. 
The intensity and length of the pulse are chosen such that the optical transition is saturated and steady 
state fluorescence is reached during the pulse. Thereby the majority of the population is transferred to 
the slowly decaying triplet sublevel, while the two other triplet levels and the ground and excited states 
are nearly empty. The laser pulse is followed by a short idle time (~50 ns) to allow any residual 
population of the excited state to decay. 

To readout the triplet spin state, the fluorescence is recorded during a laser pulse that is identical to the 
preparation pulse. Decaying fluorescence response curves are observed that encode the triplet 
populations. In particular, the fluorescence counts summed over a window at the beginning of the 
response profiles are proportional to the three triplet spin populations. This fluorescence is plotted in 
Fig. 3 and 4. In all data shown, a 200 ns integration window was used. 

Polarization and readout of the nuclear spin 

The nuclear spin is initialized by applying a long (10 u.s) laser pulse (duration longer than several triplet 
lifetimes) followed by an idle time (20 u.s) that ensures that the triplet population is fully decayed to the 
electronic ground state prior to any nuclear spin manipulation. The readout consists of a second laser 
pulse with simultaneous detection of the fluorescence, as in the case of the electron spin manipulation, 
except that now a longer (8 u.s) integration window is used to accumulate fluorescence contrast. As 
before, the integrated fluorescence encodes the spin state. 
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Fig. 1 1 Optical properties of the ST1 defect, a, Illustration of a diamond nano wire containing a single 
fluorescent point defect, b, SEM image of the [111] HPHT diamond sample with fabricated nano wires 
(diameter ~200 nm, length ~2 urn), c, confocal fluorescence microscope image of the sample. The 
marked wires contain single NV" and ST1 defects, d, fluorescence time trace recorded during application 
and removal of a non-aligned magnetic field. The fluorescence is quenched and enhanced for NV-, and 
ST1 defects, respectively, e, fluorescence emission spectrum of a single ST1 defect at room temperature 
and liquid helium temperature. A zero phonon line at 550(1) nm and a broad emission band in the red to 
NIR range is observed, f, fluorescence decay showing an excited state life time of 9.5 ns, g, second order 
photon correlation at increasing excitation power. The dip at delay t=0 reaches well below 0.5, proving 
the single center nature of the defect. At high excitation power pronounced bunching shoulders are 
observed, that indicate the presence of a long lived shelving state. 




frequency (MHz) magnetic field (mT) 

Fig. 2 1 Optically detected electron spin resonance, a, room temperature ESR spectrum at zero magnetic 
field. Resonant transitions are indicated by their zero-field splitting parameters D=1135(l) MHz and 
E=139(l) MHz. We observe an exceptionally high optical spin signal contrast of the D+E spin projection 
of about 45 %, b, Zeeman curve with an aligned axial magnetic field oriented along the defect main 
symmetry axis (35(1)° w.r.t. the [111] crystal axis); dots: experiment, lines: theoretical fit obtained from 
the triplet spin Hamiltonian H = D S z — S(S + l)/3 + E S x — S y + gS z B z , where S a are the usual 
spin operators for total Spin S = 1, g is the g -factor and B z is the axial field (see SI). The inset shows a 
schematic illustration of the crystallographic orientation of the defect. Green, blue and red arrows 
denote the [111], [110] (defect's main axis) and [001] (C 2v symmetry axis), respectively. 
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Fig. 3 1 Kinetics of the shelving state and energy level scheme, a, optical and microwave pulse sequence 
for the measurement of the shelving state lifetimes. After an optical pump pulse, a resonant microwave 
n-pulse is applied to exchange the populations of selected spin sublevels and thereby control the 
dominant decay constants, b, Decay curves of the shelving state population described by double 
exponential functions y = c — a r exp( —t/r{) — a 2 exp(— t/T 2 ). Blue: experimental data, red dashes: 
double exponential fits (see SI). For each panel, the microwave n-pulse is resonant with the denoted ESR 
transition, c, Assignment of the decay constants to the triplet energy levels. The initial populations (no 
ESR), are proportional to the decay constants. The small population of the short lived state is omitted. 
The observed decay constants imply the assignment of the energy levels denoted on the left, d, 
Summary of the triplet kinetics and electronic structure consisting of singlet ground- and excited states 
and a meta-stable triplet state. Transition rates between two states i and and the lifetime of a state i 
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Fig. 4 1 Coherent triplet spin manipulation, a, Optical and microwave pulse sequence and evolution of 
the triplet populations (here illustrated for microwaves resonant with the D+E transition), b, Rabi 
oscillations among the triplet spin sublevels. Thick blue lines show experimental data, thin red lines 
show theoretical fits with the model / = 7 — a 1 + a 2 cosD.t exp — t/r — a 3 exp(t/r), where 
t = 2(t q , _1 + r^ -1 ) -1 is a combined decay constant of the levels a and /? involved in the Rabi 
oscillations, x is the decay constant of the remaining level, a £ are fluorescence amplitudes and fi is the 
Rabi frequency. Note that for the upper and lower panel the combined decay constant is dominated by 
the fast decay of the | +) level. For the fit shown in the middle panel the last term is omitted, since the 
unaffected level is the |+) level, which (nearly) unpopulated. 
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Fig. 5 | Readout and control of a single nuclear spin in the spin-free electronic ground state, a, ODMR 
spectra of an ST1 defect with a nearby 13 C nuclear spin, under zero (upper panel) and moderate (lower 
panel) magnetic field applied along the defect's axis. E and A denote the zero-field- and hyperfine- 
splitting parameters, respectively. The hyperfine splitting is suppressed at small fields and becomes 
visible when the Zeeman energy is larger than the E parameter, b, Readout and control of the nuclear 
spin with the electronic triplet spin ancilla. The nuclear spin (blue arrows) is interrogated by optically 
pumping the electron spin (red arrows) into the triplet state and allowing the electron and nuclear spin 
to interact. The subsequent decay of the triplet state destroys the electron spin and leaves the nuclear 
spin unchanged. The bare nuclear spin is then coherently manipulated by NMR. Efficient transfer of spin 
polarization between the electron and nuclear spin in the triplet state is provided by exploiting an ONP 
as a consequence of LAC between the | -l,+l/2) and 1 0,-1/2) hyperfine levels (see text), c, Signal 
contrast of nuclear Rabi oscillations, corresponding to the product of the nuclear spin polarization and 
readout contrast, as a function of detuning from the LAC (zero detuning corresponds to the LAC). d,e, 
NMR spectrum and Rabi oscillations of the nuclear spin in the electronic ground state. For the NMR 
spectrum (d), an RF n-pulse of variable frequency is inserted between the preparation and readout laser 
pulses. The resonance dip corresponds to the Larmor precession of a bare 13 C nuclear spin. Blue: 
experimental data, Red: Lorentzian fit. For the Rabi oscillations (e) a resonant RF pulse of variable length 



is inserted between the initialization and readout laser pulses. Blue: experimental data, red: 
exponentially decaying cosine fit; decay constant 1.0(1) ms. 
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1 Introduction 

Here we present supporting material to the key results and arguments found in the main text. In particular, 
we analyze the properties of the ST1 defect that we observed through our employment of various optical 
spectroscopy and magnetic resonance techniques. From our analysis, we draw key conclusions about the 
structure of the ST1 defect and its electronic levels. These conclusions will provide the foundations for 
future studies of the formation of the ST1 defect and its properties. Our methodology draws many parallels 
with the vital early investigations of the NV center [TJ [51 [31 01 [7] • Further, we consider the mechanisms 
underlying optical nuclear spin-polarization, -readout and -decoherence in greater detail. Here we draw on 
early works on optical nuclear polarization [Hill] and 13 C relaxation in diamond [TUl [TTj . 

2 Sample fabrication and the occurrence of ST1 defects 

Arrays of diamond nano wires (length 2 m, diameter 200 nm, pitch 2 m) were fabricated into the cleaved 
[111] surface of an HPHT diamond (element six Ltd.) by reactive ion etching as described in detail in |12) . 
In brief: first, an HSQ mask is defined with e-beam lithography and a lift-off process. Thereafter, anisotropic 
etching is performed with an inductively coupled oxygen plasma. After the nano wire fabrication, the sample 
is annealed in vacuum at 850 °C for several hours and subsequently boiled in a mixture of equal amounts of 
perchloric-, sulfuric-, and nitric acid at 100 °C. 32 arrays, each containing 8080 nano wires were fabricated 
and the ST1 defects were found in all patches that were investigated. A randomly chosen area (lOOmlOOm) 
on one of the patches was intensively investigated in a home built confocal fluorescence microscope (Nikon 
LU Plan Fluor Epi lOOx, 0.9 N.A. air objective, 532 nm laser excitation) and 19 single ST1 defects were found 
within this area. This number is equivalent to approximately one ST1 defect per 800 nano-wires. Note that 
the number of NV defects in this area is considerably smaller than the number of ST1 defects (only about 5) 
and agrees roughly with the amount expected from the bulk NV- concentration. All ST1 defect fluorescence 
was observed on nano-wires. The bulk diamond underneath and the surface next to the wires was carefully 
investigated by fluorescence microscopy and no ST1 defects were observed. Note that this fact might be 
biased by the positive role of the nano-wires on the optical collection efficiency 13 . Due to the high collection 
efficiency of the nano wires, ST1 defects in the bulk would appear a factor of about 5 times dimmer than 
the defects on the nano-wires. This might render them hardly visible against the background fluorescence, 
which is about 10 kents./s on the plane sample surface and about 100 kents./s on an empty wire the present 
case. Since the nano-wires have a tendency to attract fluorescent material from the environment whose 
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emission then couples efficiently into the wires and produces background fluorescence, during the past, the 
sample surface had been cleaned several times by slowly scanning the laser over the surface (incident power 
~10 mW). This procedure bleaches most of the adsorbed molecules and removes most of the background 
fluorescence. At present we cannot exclude that this strong laser illumination in combination with ambient 
adsorbates plays a role for the occurrence of the ST1 defects. Thus, we can identify three aspects that might 
contribute to the occurence of ST1 defects. 

• The [111] crystal surface 

• The fabrication process, in particular the plasma etching using oxygen gas (trace amounts of contam- 
inats from the chamber and electrodes may also play a role) and subsequent annealing 

• Adsorbate molecules in combination with strong 532 nm laser illumination 

As we will discuss below, oxygen is a possible defect constituent since it does not have a nuclear spin and is 
capable of sp3 bonding. Future oxygen implantation and reactive ion etching studies can be used to probe 
this hypothesis. Recently a defect with a ZPL at 547.2 nm has been reported in HPHT diamonds grown 
with a Ni-containing catalyst at elevated temperatures (^1400 °C); PL spectra measured under 337.1 nm 
excitation p3]. We think that the difference of the observed ZPL compared to the ZPL observed in the 
present work at 550(1) nm is too large to identify both defects as the same species. Nevertheless, a shift of a 
ZPL e.g. due to lattice deformations or co-presence of another impurity species might still be possible. More 
importantly, for the 547.2 nm defect, an excited state lifetime of 150 s was reported which differs substantially 
from the excited state lifetime of the ST1 defect of about 10 ns. Our current working hypothesis is that the 
defect observed in the present work is a different one. 



3 Properties of the optical fluorescence 

Under 532 nm excitation, the fluorescence of the ST1 defect exhibits a zero phonon line (ZPL) at 550 nm 



(2.25 eV) and a broad vibronic band that extends to lower energy ~ 750 nm (refer to figure SI). The 
association of this fluorescence with the ST1 defect (and not with some other unidentified emission states 
in the nano-wires) was established by two different methods. Firstly, in order to verify that this emission is 
responsible for the observed single photon emission, the second order correlation measurement was repeated 
in the presence of a band pass filter centered at 550 nm (band width 25 nm) and identical anti-bunching 
behavior was observed. In order to confirm that ESR gives rise to the enhancement of this emission, the PL 
spectra were measured while applying on- and off-resonant microwave excitation. Figure [Si] shows that both 
the ZPL and its phonon side band were enhanced when the on-resonance microwave was applied. Note that 
the PL without MW was identical to PL with off- resonance MW (not shown). 

The energy of the fluorescence and the existence of the sharp ZPL implies that the electronic transition 



occurs between energetically discrete electron orbitals within the diamond bandgap (refer to figure S2| [15] . 
As the shape of the vibronic band in the PL spectra of figure [Si] is distorted by the first and second order 
Raman resonances of the 532 nm excitation, it is not possible to directly establish the Huang-Rhys factor of 
the fluorescence. The Huang-Rhys factor describes the relative intensity of the ZPL and the vibronic band 
and is a measure of the lattice relaxation associated with the electronic transition |15j . Since the Raman 
resonances should not change in the presence of microwaves, the difference in the PL spectra while applying 
on- and off-resonance microwaves provides a reasonable picture of the ZPL and vibronic band that is free of 
distortions from the Raman resonances. It is apparent that like many other color centers in diamond, the 
Huang-Rhys factor at room temperature is appreciable (estimate at least greater than one) and that there is 
an appreciable lattice relaxation associated with the electronic transition p] [16] . The difference between the 
first moment of the vibronic band and the ZPL energy is ~ 0.15 eV at room temperature, which provides an 
estimate of the Anti-Stokes shift and lattice relaxation energy associated with the transition [15 . Further 
studies of the vibronic band will provide additional information on the electron-phonon interactions of the 
ST1 defect Q7]. 
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Figure SI: Enhancement of emission induced by ESR. Low panel: Two PL spectra measured at room 
temperature with on- and off-resonance continuous ESR excitation. Upper panel: Difference between two 
spectra. Enhancement of both ZPL and phonon side band can be seen. The sharp peak at 573 nm is the 
Raman peak. 
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Figure S2: Electron orbitals of the ST1 defect. The shaded regions represent the diamond valence 
and conduction bands separated by the indirect bandgap of ~ 5.5 eV. The solid lines within the band gap 
represent the energies of the orbitals of the ST1 defect that are involved in the optical transition. The 
a orbital is the highest occupied orbital in the electronic ground state. Optical excitation results in the 
promotion of an electron from the a orbital to the lowest unoccupied orbital b. The symmetries of the 
orbitals are unknown. There may exist other orbitals of the ST1 defect in the bandgap. 
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Figure S3: Rate equation model. The ground and excited states are denoted G and E. The three triplet 
spin levels are denoted |0), |— ) and |+). Rates and lifetimes are denoted by greek letters 7 and r, respectively. 

4 Electronic structure and optical cycle 

It was discussed in the main text that the key experimental observations that suggest the existence of a 
metastable triplet state are the increase of the fluorescence under application of a strong misaligned magnetic 
field or resonant microwave field, and pronounced photon bunching under saturated optical excitation (refer 
to Fig. 1, main text). Such behaviour is well known for fluorescent molecules with singlet ground and 
excited states and a long lived excited triplet state that features different lifetimes of the individual triplet 
spin sublevels [IS]. The triplet state plays the role of a shelving state. When the population is cycled 
between the ground and excited state by optical pumping, part of the population gets trapped in the triplet 
state. This Inter-System-Crossing (ISC) leads to an overall reduced fluorescence since the ISC transitions are 
typically non-radiative and the lifetime of the triplet state is long (ISC transitions are only weakly allowed 
by spin-orbit or spin-spin interaction). The long lifetime leads to characteristic 'bunching shoulders' in the 
second order photon correlation. When either the population- or the decay rates of the triplet sublevels 
are different, the sublevels are not equally populated, i.e., the spin becomes polarized. Spin polarization 
is typically built up during several optical pump cycles. When the decay rates of the triplet sublevels are 
different, a mixing of the triplet populations (either by a resonant microwave field or by a strong misaligned 
magnetic field), results in a faster or slower deshelving of the triplet state. This leads to an increase or 
decrease of the emission intensity. 

In the following, we analyze the optical cycle and provide conclusions about the defect electronic structure. 
We first present a 5-level rate equation model to describe the optical cycle. Thereafer we experimentally 
determine the transition rates. Finally we draw conclusions for the electronic structure and symmetry of the 
electronic states. 

4.1 Rate equation model 

The rate equation model is illustated schematically in figure |S3| It consists of five levels: the ground and 
excited states G and E and the three triplet spin levels |0), |— ) and |+). Rates between levels and v 
are denoted by 7^ and the lifetime of a levels is denoted by r. 7gb is the optical pump rate, jeg is the 
spontaneous emission rate, Jeu, v = |0), |— ), |+) are the triplet population rates and 7c = |0), |— ), |+) 
are the triplet decay rates. The rates are related to the lifetimes r = lv The time evolution of the 
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populations n — (nc, tie, «|+) , n \-)i n |o)) i s given by the rate equations 



nG=- IGETIG + lEGflE + l\+)G n \+) + 7|->G«|-) + 7|0>G«|0) 

n E = + iGEna - Ieg^e - (7e|+) + 7e|-> + 7b|o))"b 
™|+> = + 7b|+)«b ~ 7|+)G«|+) (S.l) 
«|-> = + 7e|->«£ - 7|->g«|-> 

»|0) = + 7b|o)?^b - 7|o)G"|o> 

Using the conservation of the total population 

1 = tig + n.E + n\+) +«|-> + n |o)j ( s - 2 ) 



one euqation and one variable in Eq.(S.l) can be eliminated. The result is a regular first order linear 



differential equation with constant coefficient matrix A and right hand side vector b. 

h = An + b, (S.3) 

which has the solution 

n(t) = e At n + A- 1 e At b- A'% (S.4) 

where no are the initial populations. Using such a model, all photophysical properties, such as the satura- 
tion behavior, anti-bunching, fluorescence response to an optical pump pulse, etc., can be easily computed 
numerically by diagonalizing the coefficient matrix A or in some cases analytically. In the following sec- 
tions, we will experimentally determine the decay rates 7|+)g,7|->G)7|o>G an d the total population rate 
Iet = 7.E|+) + 1e\-) + 7b|o) °f the triplet level, where we assume that the population rates of the triplet 
state are faster than the decay rates and are equal for all spin sublevels, and the triplet decay rates are 
distinct. In this case, the populations of the triplet levels under saturated optical excitation are proportional 
to their lifetimes. 

4.2 Triplet lifetimes 

As discussed in the main text, in order to determine the triplet decay rates, the recovery of the ground 



state population was observed by a pulsed optical experiment (refer to Fig. 3, main text and figure S4|. 
First a long laser pulse is applied (duration longer than the triplet lifetimes, power high enough to saturate 
the optical transition). At the end of the pulse the system is in the steady state under saturated optical 
illumination, i.e., the triplet levels are maximally populated. The relative populations of the sublevels are 
proportional to their lifetimes. Next, the triplet levels are allowed to decay for a variable amount of time in 
the range few ns to several s. Finally, the population of the ground state is measured by applying a second 
laser pulse that is identical to the first pulse and simultaneously recording the fluorescence. The integrated 
fluorescence is proportional to the ground state population before the pulse. Most data was accumulated 
until the standard deviation of the measured ground state populations (given by photon shot noise) was 
below 1%. Exponentially increasing fluorescence curves were obtained. The time constants are given by 
the lifetimes of the triplet sublevels. The data was fit with single-, double- and triple-exponential models 
and the goodness-of-fit Q — r(0.5^, 0.5xo) was evaluated for each fit [19]. Here T is the complementary 
incomplete gamma function, v is the degrees of freedom and Xo is the x 2 evaluated at the minimum. Single 
exponential fits typically result in goodness of fit values Q < 1CU 10 and can be discarded. Most Double- and 
triple-exponential fits result in goodness-of-fit values in the range 1 < Q < 1CU 3 and are both acceptable. 
The triple-exponential models typically contain one component whose amplitude is much smaller than the 
two other components and therefore this component can be safely neglected. This reflects the fact that 
acceptable goodness-of-fit values are readily obtained with double exponential fits. We conclude that double 
exponential fits are faithful models for all data. To assign the lifetimes to the spin sublevels, we perform the 
measurement with an additional microwave 7r-pulse inserted after the optical pump pulse. The microwave 
power is high enough, such that the duration of the microwave pulse is short compaired to the lifetimes. The 



18 



a 



init 



count 




90 
80 
70 
60 
65 

60 
55 
50 
45 
40 
35 
30 



1 1 1 1 1 1 1 

4ji ii|i>i'li>ii "HV t » M ii 11' i ri" i ' , r't i ' i i r , 'f Ttt i nr'rrt- 



no MW 



D + E 



D-E 





— i — i — i 


— i — i — i — 






2E 




i i i 





6 8 10 12 14 16 

time (|js) 



Figure S4: Experimental determination of the Triplet lifetimes, a, Pulse sequence, b, Measured 
decay curves for defect C (see text), without microwave pulse and with a 7r-pulse tuned in resonance with 
the three ESR transitions. Blue dots: experimental data, solid green lines: single exponential fits, solid red 
lines: double exponential fits. 
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Table SI: Measured triplet lifetimes of three STl defects A,B,C. The data in each cell shows (from 
top to bottom) the two time constants [ns], the two amplitudes [counts] and the goodness-of-fit of a double 
exponential model. The last row summarizes the observed time constants [ns]. The mean and standard 
deviation over equivalent time constants yields r|o) = 2500(204) ns, T]_) = 1072(350) ns, t]+) = 209(44) ns 



microwave frequency is tuned in resonance with each of the triplet spin transitions. The microwave 7r-pulse 
swaps the population of the triplet sublevels and thereby exchanges the amplitudes of the corresponding 
exponential decay curves. Table [ST] summarizes the measured decay constants and amplitudes and specifies 
the goodness-of-fit values for three different STl defects. Importantly, there are only minor deviations among 
the different defects. Generally three distinct time constants are observed, of about 200 ns, 1000 ns and 
2500 ns. The result is summarized in the last row. The result is used to assign the time constants to the 
spin levels as shown in the illustration below the table. To obtain an overall average for the decay constants, 
we evaluate the mean and standard deviation over all equivalent time constants, yielding T| ) — 2500(204) 
ns, T|_) = 1072(350) ns, T| +> = 209(44) ns. 



4.3 Triplet population rates 



In the following, we determine the triplet population rates by two different methods. One approach exploits 
the second order correlations, </ 2 \ of the photon emission at different excitation powers. The numerical 
solutions is fit based on the five level rate model described in equation (S.ll yiel ding the population rate 
parameters. The measured g^ 2 ' functions and fitting results are shown in figure S5 Two representative 
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Figure S5: Determination of the population Rates. Two upper panel: measured functions at 
low and high optical excitation powers. Red lines are the fitted curves. Instead of the excitation powers, 
the extracted pumping rates are shown. Lower panel: The obtained population rates as a function of the 
pumping rate. Only a common value 5.9(6) MHz was obtained because the three population rates could not 
be distinguished. Blue dotted line indicates this value which is almost constant. 
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time offset [ns] 

Figure S6: Determination of the overall triplet population rate by pulsed optical excitation. 

a) Schematic of the measurement and data analysis. The stimulation consits of optical pulses that are 
separated by idle times. Both are longer than the triplet lifetimes. The resulting fluorescence response 
profile is processed by integrating the fluorescence over a window of fixed length and moving this window 
over the profile, starting from the initial rising edge, b) Experimental data (blue solid line) and double 
exponential fit (red dashed line), providing two time constants t\ = 85ns and T2 = 724ns. 



g( 2 ' functions measured at low and high excitations are shown. Both the anti-bunching part and bunching 
shoulder can be well identified. Since the excited singlet state life time and the triplet state life times 
are already known, the measured </ 2 ' functions can be fitted with four free parameters, namely the pump 
rate and the three population rates. Three population rates, however, could not be distinguished within 
the error range, instead the fit was simplified to two free parameters by considering the total rate jet = 
(56(7) ns) -1 , where we assume Je\+) = TBI—) = Tbio)- Note that whilst ISC selection rules likely imply that 
the populations rates could differ, such differences will not significantly effect the application of the rate 
model in the present section. 

A second method for determining the triplet population rates is to use a pulsed optical experiment. We 
expect that the population rates of the triplet sublevels are similar for all three levels and the time scales 
of population and decay are well seperated (the population rates and decay rates differ by about an order 
of magnitude) . If this was not the case, we would expect a more complicated decay behavior in the lifetime 
measurements. In this case, the total transition rate from the excited state to the triplet levels can be 
determined from the fluorescence response of the system to a single optical pulse as shown in Fig S6 1. The 
length of the laser pulses and the idle time between two pulses is much longer than the triplet lifetime. If the 
laser is set to high power such that the optical transitions is saturated, the ground state is essentially depleted 
and first all population resides in the excited state. The population then starts to leak to the triplet state. 
Thus, the length of the excitation pulse determines how much population is transferred. Since the lifetimes 
of the triplet state are much longer than the population rates, recycling of population through the triplet and 
ground state during this process is negligable. Therefore, an approximate exponential increase of the triplet 
population is expected with time constant tet = Iet- As before, the triplet population can be measured 
by a readout laser pulse and recording the fluorescence response profile. In the present measurement, the 
initialization and readout laser pulses can be merged to a single pulse and the distinction between the two 
pulses is achieved by moving the integration window for the determination of the triplet population over 
the pulse profile (refer to figure S61). The result of this measurement is shown in Figure S6d. A double 
exponential decay is observed. The short time constant reflects the buildup of the triplet population. The 
longer time constant reflects the time required to reach the steady state under optical pumping and is given 
by a combination of the different decay rates. The combined population rate is found from the short time 
constant as jet ~ (85 ns) -1 . This value agrees reasonably well with the value of about 56 ns obtained from 
the photon auto-correlation. Appreciating both independent methods, we estimate jet ~ (70(20) ns) -1 . 
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4.4 Application of the rate equation model 

In this section we provide further detail about the rate equation model. In particular, we present that the 
recovery of the ground state population as observed in Fig 3 in the main text is the evidence for the triplet 
meta stable state and analyze how the observed amplitudes of exponentially decaying components of the 
triplet levels are related to the rate parameters. To this end, the analytical solutions for the ground state 
spin recovery without and with 7r-pulses will be introduced. The steady state solutions based on the rate 
model can be easily found (not shown) and can be considered as initial values at t = for the ground 
state recovery measurement. Then the spin populations of each states after time r' (while jge — 0) can be 
found. In the experiments, a 6 s long readout pulse was applied at t = r' and the fluorescence response was 
integrated only over a short (20 <~ 200 ns) period at the very beginning of the readout window. Because no 
significant difference in the integrated signal was found up to a 200 ns long integration window, we suppose 
that the repolarization hardly affected so that the integrated intensity is directly proportional to the ground 
state spin population after r'. Thus the ground state spin population at t = t' is only of our interest, which 
is 

n G {r') = & ■ n G (Q) ■ e" rV - V — — • m(0) ■ e^^' + C • n G (0) (S.5) 

H+>.l->.|0) "° 

where rij(0) is the initial value for each state which is equivalent to the steady state value, T — Jeg + 1e\+) + 

1e\-) +7b|o>, - + 2-,i=|+>,|->,|o> r*.(r- 7iG ) ' ancl °° ~ 1 Cl + ^=|+>,|->,|o) liG -(~, EG - 7iG ) • 

When an ideal resonant 7r-pulsc between two sublevels, |+) and |— ) for instance, is applied right after the 
initialization pulse, the initial population of these two states are exchanged. The solution for this case 
becomes 

n G ,\ + )\-)(r') = (7 1 .n G (0)-e- rV - £ K(0) + ■ n,(0)) • e~^' 

-7^— •n|o)(0)-e-^ GV +C o ,| +) |_ ) -n G (0) (S.6) 

1 _ 7|0)G 

where Co,|+>| > i s a constant similar with Co- Solutions for the two other cases nG,|+)|o>( 7 "') an( i n G,|-)|o) ( T ') 

also can be found (not shown). Both solutions consist of four exponential terms, but the first term has little 
contribution because its decay is faster than the others (r > jeg S> 7jg) by at least one order of magnitude, 
and its amplitude is very small too as will be seen later. In addition, because the found triplet state life times 
are at least one order of magnitude slower than the excited singlet state life time, we get 7j G /(r — 7, G ) *C 1 
and r/(r — 7ic) ~ 1- Thus the above solutions can be rewritten as 

n G (r')^- Yl rn(0) ■ e~™- T ' +C a -n G (0), (S.7) 

i=|+>,|->,|0) 

n G ,| +> |_>(/)»- Yl n,(0)- e -^ V -n| O) (0)- e -^>« V +C o ,| +) |_ ) ■ n G (0). (S.8) 

»>j'=I+>,|->,*^j 

Thus one can conclude that the observed ground state recovery curve can be explained well by three 
exponential decay terms of which decay constants are the decay constants of the triplet sublevels or even the 
double exponential decay if one of the three decay terms is too small. And the inserted 7r-pulse does not alter 
the decay constants but only exchanges the amplitudes of the sublevels under the influence of the applied 
7r-pulse as long as the excited singlet state life time is faster than the triplet state life times. Therefore, one 
can assign the measured decay constants to each sublevels by comparing the amplitudes of each decay terms 
obtained with and without applying 7r-pulses. In order to make an indisputable conclusion, it is better to 
compare the numerically simulated amplitudes of each decaying components with the observed ones. 
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Figure S7: Numerical simulations of the ground state recovery measurements. Red lines are the 
simulated recovery curves, cyan, blue, and green curves are the exponential decays of |+), |— ), and |0), 
respectively. The first exponential terms (violet lines) in equation (S.5| and (S.6) are almost ignorable 
because of their fast decay and small amplitudes. 



The numerical simulations using the exact solution, equations (S.5) and (S.6 1 with the set of parameters 
of 7-EG = (10ns) _1 
(1000 ns)- 1 



7|0)G 



7gb = (10 ns) \ 7 B | + ) = 7£|_ 
(2500ns) -1 are shown in figure 



S7 



7b|o> = (170ns)- r 7^ +)G = (200ns)- 1 , 7 |_ >G = 
Note that the common value of 170 ns was used 



for three population rates and the typical value for the optical pumping rate jge = (10ns) -1 was used. 
The first term in equations (S.5) and (S.6) (violet curves in figure S7) of which decay constant is T = 
Ieg +7b|+> +7b|-) +7b|o) ; ls negligible because it decays very quickly and the amplitude is minute. Among 
the remaining three exponential decays originated from the decays of the triplet sublevels, the amplitude of 
the fastest decay is smaller than others, so this decay is hardly detectable in measured data as already shown 
above. Overall, these simulated curves are in close agreement with the measured ones and can explain the 
measured curves very well. The simulated and measured amplitudes of each decay terms are summarized 
in table |S2| The measured amplitudes without applying 7r-pulses are in very good agreement with the 
simulations. The exchange of the initial populations by the coherent 7r-pulses are well described in both 
measured and simulated values, and they are in good agreement. The found discrepancy between the size 
of the simulated and the measured amplitudes can have a few origins: intersystem crossing among each 
sublevels, non-ideal 7r-pulses. 
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Table S2: Comparison of the measured amplitudes with the simulations. The observed switching 
of the initial populations by 7r-pulses are in consistent with the simulations. 



4.5 Conclusions for the electronic structure 

As discussed in the previous sections, the observed behaviour of the optical fluorescence implies that the 
optically detected magnetic resonance occurs within a metastable spin triplet level that exists energetically 
between the ground and optically excited levels of the defect. It was also discussed in the main text that 
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Figure S8: Observed electronic structure of the ST1 defect. Levels are labelled by their appropriate 
term symbols. Note that the orbital symmetries of 1 E and 3 T are unknown. The vertical arrow denotes the 
optical ZPL at 550 nm (2.25 eV). The estimate ~ 0.15 eV obtained for the lattice relaxation energy in the 
previous section can be combined with the ZPL to estimate an energy of ~ 2.40 eV between the 1 Ai and 
1 E electronic levels. Other electronic levels may be present, but are yet to be detected. 



the nuclear resonance of the 13 C nuclear spin performed whilst the defect is in the electronic ground state 
implies that the ground electronic level is a spin singlet. These facts can be used to conclude various aspects 
of the defect's electronic structure. 

According to Hund's rules, a spin singlet ground state can only occur if the ground electronic configuration 
is closed (i.e. filled electron orbitals) [20]. Using the electron orbital notation of figure S2 the ground 
electronic configuration is a 2 , where all lower energy orbitals (i.e. orbitals of the valence band) are taken to 
be filled. Applying Unsold's theorem, a closed configuration is necessary orbitally symmetric (denoted Ax) 
[20] . such that the term symbol of the ground electronic level is 1 A%. Optical excitation promotes one of 
the electrons to the higher energy orbital 6, such that the configuration of the optically excited level is ab. 
Spin selection rules of optical transitions imply that the optically excited state is also a spin singlet. As the 
symmetry of the orbitals are unknown, it is not possible to determine the orbital symmetry of the optically 
excited state, so for convenience it will be denoted by the term symbol 1 E. The spin triplet level may be of 
the same configuration as E (the triplet being lower in energy in accordance with Hund's rules) or another 
excited configuration. Consequently, the orbital symmetry of the triplet level is also unknown and it will be 
labelled by the convenient term symbol 3 T. There may also exist other electronic levels between the triplet 
and singlet levels and also above 1 E (which are relevant as they may be populated under 532 nm excitation). 
However, such additional levels have not been detected thus far. The aspects of the electronic structure that 
can be concluded from our observations are depicted in figure [S8| Future optical polarization and stress 
studies can be used to determine the orbital symmetry of the optically excited singlet level [2]. 

Given the observed aspects of the defect's electronic structure, a working model can be established of 
the defect's optical cycle (refer to figure S9). Beginning in the ground state, optical excitation transfers the 
defect into a vibronic level of E, which is followed by rapid non-radiative decay through the 1 E vibronic 
levels, such that thermal equilibrium is established within the level's lifetime. From the X E level, the defect 
may decay radiatively back to the ground state or non-radiatively via an intersystem crossing (ISC) to the 
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Figure S9: Model of the optical cycle of the ST1 defect. Known electronic levels are denoted by 
solid lines. The continuum of vibronic levels that give rise to the optical absorption and emission bands 
are depicted as shaded regions. Solid vertical arrows denote radiative transitions. Dashed arrows denote 
non-radiative decay within the vibronic levels of 1 A% and 1 E as well as the ISCs between the singlet and 
triplet levels. 



metastable triplet level. The defect decays from the 3 T level by a subsequent ISC to the ground state. The 
ISCs are allowed via spin-orbit coupling of the triplet and singlet levels and are mediated by phonon emission 
[3]. The nature of spin-orbit coupling implies that the ISCs will be to some degree dependent on triplet spin- 
projection. The net effect being different storage times in the triplet level for different spin-projections that 
results in a modulation of the optical fluorescence if the defect is transferred between the spin-projections 
of the triplet. The presence of other electronic levels in the optical cycle will to some degree change the 
behavior of the cycle from that described by the model depicted in figure |S9| No such variance has been 
detected to date. 

5 Triplet level spin-Hamiltonian, zero-field splittings and spin axes 

The spin-Hamiltonian that is appropriate to model the fine structure of the triplet level takes the general 
form El] 

H = S-T>- S + S-g- B (S.9) 

where S is the dimensionless 5 = 1 spin operator, B is the magnetic field at the defect, and D and g are 
each 3x3 tensors to be determined experimentally. The first term describes the effects of the crystal field 
that typically arise from electronic spin-spin interactions and/ or second-order spin-orbit interactions. The 
second term describes the electronic Zeeman interaction as governed by the tensor g, which typically differs 
from the isotropic g-factor of a free-electron g e = 2.0023 by orbital angular momentum and spin-orbit effects. 

By an appropriate choice of the defect coordinate system (xyz), the crystal field tensor D may be 
diagonalized. Furthermore, as our experiments did not detect any anisotropy in the electronic Zeeman 
interaction, the tensor g may be reduced to an isotropic g-factor. The simplified spin-Hamiltonian is thus [3J 

03] 

H = D[S 2 Z - S{S + l)/3] + E{S 2 X - S 2 ) +gS-B (S.10) 
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Figure S10: Triplet level zero field fine structure. Note that D and E have been assumed positive. In 
terms of the spin states \m a ) with triplet spin-projection m s , the zero field eigenstates arc |+) = (| + 1) + 
|-l»/v / 2,|-) = (| + l)-|-l»/v / 2and|0). 
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Figure Sll: Distributions of the ZFS of 19 ST1 defects. The variance of the parameters is consistent 
with a distribution of strain within the nano-wires. 



where D and E are zero field splitting (ZFS) parameters and g = 2.0(1) is the observed isotropic g-factor. 
The zero field fine structure of the triplet level that is described by the above spin-Hamiltonian is depicted 



in figure S10 The distributions of the observed ZFS parameters of 19 defects are shown in figure Sll 



The relatively small variance of the parameters strongly indicates that they are intrinsic properties of the 
ST1 defect that are slightly modified by the local strain of each defect site. Note that the imprecision of 
the observed g-factor is due to the experimental uncertainty of the magnetic field calibration. Also note 
that a precise measurement of the g-factor could give additional insight into the symmetry of the defect: 
Anisotropy of the g-factor arises due to spin-orbit coupling between the electronic levels and these couplings 
obey symmetry selection rules. C2v symmetry implies negligible spin-orbit coupling and thus likely an 
isotropic g-factor, while in case of a symmetry lower than C2v, an anisotropic g-factor could be present. The 
precise sign of the ZFS parameters can not be determined from our measurements. 

The spin transition frequencies for different orientations of a static magnetic field are depicted in figure 



S12 The spin resonances have clear relationships when the magnetic field is aligned along one of the defect's 
coordinate axes. The orientation of the magnetic field can be simultaneously measured with respect to the 
crystal coordinate system using three inequivalent nearby NV _ centers and knowledge of their [111] symmetry 
axis (method described elsewhere [3T]). Consequently, the ST1 defect's major (z) axis was determined by 



27 




a 

1400 
1200 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

{rad) 

Figure S12: Magnetic field alignment and major symmetry axis. The plots depict the variation of 
the three spin resonances as a function of magnetic field orientation: (a) B x — Bs'md, B y — 0, B z — Bcosd; 
and (b) B x = Bcos<p, B y = Bsin^, B z = 0. The values D = 1128 MHz, E = 139 MHz and B = 300 MHz 
were used to produce the plots. Relationships between the resonances are evident when the magnetic field 
is aligned with the defect's z axis (9 — 0), x axis (<fi — 0) and y axis (</> = tt/2). (c) Major symmetry axis of 
the ST1 defect. Blue lines indicate the four NV axes in the diamond crystal. 



monitoring the ODMR frequencies of a chosen ST1 defect and nearby NV centers whilst manipulating 
the orientation of the magnetic field generated by an electromagnet placed above the sample. The elicited 



major axis together with the four NV~ orientations is shown in figure S12:. The major axis is on a plane 
defined by two NV~ axes (NV1 and NV2 axes in this plot) and off by 35(1) ° with respect to these two axes 
and perpendicular to the other two NV - axes. These findings lead to the conclusion that the major axis is 
approximately in the [110] crystal direction. The same conclusion has been obtained from several other ST1 
defects. 

The ZFS parameters of molecule-like defects in diamond, such as the NV - center, are predominately 
due to spin-spin interactions, whereas the ZFS parameters of atom-like defects, such as transition metal ion 
impurities, arise largely from spin-orbit interactions [5J[inj- Assuming that the ZFS parameters arise from 
spin-spin interaction between the two unpaired electrons of the triplet level, the ZFS parameters provide 
insight into the distribution of unpaired spin density associated with the triplet level via the following 
expressions [15] 

3„ 22 o ,l-3(*? a /r? 2 ). 



D =2^4. 



12 

,,2 



where o is the vacuum permeability, b is the Bohr magneton, X12, yi2, z\2 are the coordinates of the 
displacement vector r\i with magnitude r\i that connects the positions of the two electrons, and the angle 
brackets denote the direct integral over the unpaired spin density. Note that the exchange integral is neglected 
in the above. 

The interpretation of the observed ZFS parameters of the ST1 defect using the above expressions indicates 
that the unpaired spin density is elongated along the major (z) axis of the defect and is not symmetrical 
in the x — y plane of the defect |15j . To explain this interpretation, an analogy may be drawn with the 
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NV~ center, where the ground state unpaired spin density is equally distributed over the three dangling sp 3 
bonds of the carbon atoms surrounding the vacancy [3] . The trigonal symmetry of the NV~ center implies 
that its ground state spin density is symmetrical in the (111) plane in the absence of strain, but is elongated 
along the [111] direction. The E parameter of NV~ is therefore zero in the absence of strain, whilst the D 
parameter is non-zero. In the presence of strain, the trigonal symmetry of the NV _ center is lowered and its 
ground state spin density is distorted in the (111) plane, thereby resulting in a non-zero E parameter. The 
non-zero E parameter of the ST1 defect implies that it is a defect of less than trigonal symmetry [51 115j. 
Combining this conclusion with the orientation of the defect's major axis along the [110] direction restricts 
the set of possible point group symmetries of the ST1 defect to just including C<i v and its sub-groups: C2, 
C lh , and C 1 . 

The value of D for the ST1 defect may also be used to crudely estimate an upper bound of 4 angstroms 
for the average separation (ria) of the two unpaired electrons [T5]. A similar estimate applied to the 
ground state of the NV~ center yields 3 angstroms [5]. By comparison, it can be concluded that the spin 
density of the ST1 defect is likely to be distributed over a greater region than the ground state spin density 
of the NV _ center. The ZFS parameters are thus consistent with an unpaired spin density distributed over 
the orbitals of several atomic sites along the [110] crystal direction. The ZFS parameters are thus consistent 
with some degree of preferential distribution of the unpaired spin density over the orbitals of several atomic 
sites along the [110] crystal direction. 



6 Triplet hyperfine structure 

Two of the ST1 defects studied exhibited hyperfine structure of the triplet level due to magnetic interactions 
with a 13 C nucleus. One of them showed a particularly large hyperfine splitting greater than 100 MHz and 



was used throughout this work (refer to figure 4 main text). Figure S13 shows the hyperfine splittings of 



this defect in the vicinity of a level-anti-crossing between the |0) and | — 1) electronic spin manifolds. These 



interactions are described by the addition of the following term to the triplet spin-Hamiltonian (S.10) [3, 15 



V hf = S-A-I (S.12) 

where / is the dimensionless 1=1/2 nuclear spin operator and A is a 3 x 3 tensor to be determined 
experimentally. With an appropriate choice of nuclear spin coordinate system ( XY Z) , the tensor A may be 
diagonalized, simplifying the above to [15] 

V h f = A XX S.J X + AyySylY + A ZZ S Z I Z (S.13) 

Note that the nuclear spin coordinate system (XYZ) is in general different to the defect's coordinate system 



(xyz). Figure S14 depicts the triplet hyperfine structure at zero field and in the presence of a static magnetic 
field aligned along the defect's major axis (B z ) (nuclear Zeeman interaction ignored). At zero field, the 
hyperfine structure is suppressed by the effects of the crystal field parameter E (refer to figure 4 main text). 



The observed hyperfine resonances as functions of B z (data shown in figure S13), together with their least 
squares best fit obtained using the spin-Hamiltonian model, are depicted in figure |S15| The fitted values of 
the various spin-Hamiltonian parameters are contained in table |S3| Note that no difference could be detected 
between the hyperfine parameters A xx and A yy , such that they both can be equated to a single parameter 
A , - A -A 

^xx — -tiyy ■ 



D (MHz) 


E (MHz) 


A zz (MHz) 


A± (MHz) 


1134.7(7) 


139(2) 


-117(1) 


-94(2) 



Table S3: Triplet level hyperfine parameters. The parameters are those of the single ST1 defect whose 
hyperfine structure was studied in figure S15 Note that A xx = A yy = Aj_. 



The hyperfine parameters can be used to interpret the character of the unpaired electronic spin density 
distribution at the 13 C nucleus. The magnetic hyperfine interaction contains distinct Fermi contact / and 
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Figure S13: Hyperfine structure at the LAC. Upper panel: D + E transition, lower panel: D — E 
transition. The color maps show experimental ODMR data. White lines show theoretical results. The 
approximate eigenstates linking the transitions arc denoted on the sides of the plots. 
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Figure S14: Triplet level hyperfine structure. The zero field energies are correct to first order in the 
hyperfine interaction. The approximate eigenstates in the presence of an aligned magnetic field fulfilling 
B z 3> E are expressed in terms of the electronic spin states (|m s )), where m a denotes electronic spin- 
projection, and the nuclear spin states (| t)? 
spin projections. 
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dipolar d terms, such that [3J [5] 



A zz = f + d 

A± = f-2d (S.14) 



Using the fitted values of A zz and A± contained in table S3 the values / = — fOf (2) MHz and d = —7.8(9) 
MHz can be determined. The electronic orbitals ip that comprise the unpaired spin density associated with 
the 13 C atom can be well described by linear combinations of 2s (</> s ) and 2p (<f) p ) atomic orbitals [3j [5] 

ip = c s (f> s + c P 4>p (S.15) 

where c s and c p are coefficients that satisfy the normalization condition f = |c s | 2 + |c p | 2 . The Fermi contact 
and dipolar terms may be expressed in terms of the orbital coefficients and the portion of unpaired spin 
density rj contributed by the atomic orbitals of the 13 C atom via 3, 5 J 

Sir 

f = y ^ffeSff™n|c s | 2 ^|0 s (O)| 2 



2 1 

d= ^^9eBgnn\Cp\ 2 r]((f) p \-^\(f) p ) (S.16) 



where g n is the 13 C nuclear g-factor, „ is the nuclear magneton, </> s (0) is the value of the 2s orbital at 
the nucleus and r is the distance of the electron spin from the nucleus. As per established practice in the 
electron spin resonance of defects in diamond [3[ Q21 G3 [55], the values of |^ s (0)| 2 and (</> P |^|0p) obtained 
from self-consistent field calculations of atomic carbon can be used in conjunction with the observed values 
of / and d to estimate |c s | 2 , |c p | 2 and r\. The estimated values are 0.27, 0.73 and 0.1, respectively. 

The estimated values of |c s | 2 and |c p | 2 are very close to the values of a standard sp 3 bonding orbital of 
diamond, being 0.25 and 0.75, respectively. Consequently, it appears that the spin density at the 13 C site 
of this particular ST1 defect occupies a typical sp 3 orbital. Furthermore, only ~ 10% of the spin density is 
associated with that site, implying that the spin density is distributed over a number of atomic sites. This 
conclusion, when combined with the observed ~ 1 MHz ODMR linewidths and typical hyperfine interaction 
energies, implies that if the defect possessed an intrinsic nuclear spin it would most likely have been detected. 
Hence, it may also be conservatively concluded that the defect does not possess an intrinsic nuclear spin. 

7 Defect structure 

Summarising previous sections, the properties of the ST1 defect structure that we may conclude from our 
observations are: 

• The optical transition occurs between energetically discrete electron orbitals within the diamond 
bandgap and is accompanied by an appreciable lattice relaxation. 

• The ground and optical excited electronic levels are spin singlets 

• There exists a metastable spin triplet level between the singlet levels 

• The triplet level's major axis is in the [110] direction. 

• The defect has C2 V or lower point group symmetry. 

• The defect does not contain an intrinsic impurity that possesses a nuclear spin. 

• The mean separation of the unpaired electrons of the triplet level is less than 4 angstroms. 

• The unpaired spin density is distributed over a number of atomic sites. 
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Figure S16: Possible regions of the diamond lattice occupied by the ST1 defect. The diagrams 
depict two portions of the diamond crystal structure: (a) a portion centered on an atomic site and (b) a 
portion centered on the midpoint between two atomic sites along the [110] direction. The transparent sphere 
in each diagram has a diameter of 4 angstroms and indicates the region where the unpaired spin density of 
the triplet level is localized. The red arrow denotes the C2 symmetry axis of the [110] direction and the blue 
arrow is parallel to the [110] direction. The plane formed by the C2 symmetry axis and the [110] direction 
is the reflection plane of the diamond lattice corresponding to the [110] direction. Black spheres represent 
atomic sites and white tubes denote sp 3 bonds. 



• 10% of the unpaired spin density is located at one particular carbon site and has the character of a 
sp 3 orbital. 



Figure S16 depicts two portions of the diamond crystal structure: one portion is centered on an atomic 
site and the other portion is centered on the midpoint between two atomic sites along the [110] direction. 
The transparent sphere in each diagram has a diameter of 4 angstroms and thus indicates the region where 
the unpaired spin density of the triplet level is localized, if the defect is atom centered or midpoint centered. 
Consistent with our observations, it is clear that for both cases, there exist several atomic sites (up to 7) 
within the region where the spin density is localized. 

The conclusion that the defect does not contain an intrinsic impurity that possesses a nuclear spin restricts 
the set of possible constituents of the defect. The restricted set contains substitutional and interstitial 
impurities without nuclear spins, as well as the structural defects of diamond: the vacancy and self-interstitial. 
The defect is likely to be a complex of these possible constituents, although given the relatively small region 
where the spin density is localized, it is probable that it contains only up two or three. The simplest structure 
of the defect that is consistent with our observations is a pair of constituents occupying atomic or interstitial 
sites along the [110] direction. A pair of the same constituent will form a defect with the highest possible 
C% v symmetry if they are both occupying atomic or interstitial sites. The pair of different constituents will 
form a defect with C'ih symmetry. More complicated complexes can form defects with C2 and C\ symmetry. 

Due to the involved fabrication process of the nano-wires, it is not possible to further restrict the set 
of possible constituents of the defect. Systematic formation studies are required to gain more information 
about the constituents. Given that oxygen plasma was used during the reactive ion etching process, it is 
capable of sp 3 bonding and the abundant le O isotope does not have a nuclear spin, the presence of le O in 
the defect is certainly possible and warrants further investigation. The significant strain that is likely to be 
present in the nano-wires potentially plays a significant role in the formation energetics of the ST1 defect 
and may explain why the defect has not been previously observed in bulk diamond. 
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Figure S17: The triplet level hyperflne level anti-crossing. Solid lines denote the hyperfine energy 
levels in the region of the LAC. The approximate states and electro-nuclear spin projection m of the levels 
before and after the LAC are depicted in black. The mixed states of the m = —1/2 levels at the LAC 
are depicted in blue. The dashed circle highlights the minor LAC between the m = +1/2 and —3/2 levels 
discussed in the text. 



8 Nuclear spin polarization and readout 

The hyperfine structure of the ST1 defect discussed in section [6] exhibits a level anti-crossing (LAC) in 
the region of a static magnetic field of strength B z ~ 42 mT that is aligned along the defect's major axis 



(refer to figure S17). Since B z and D are the dominate terms of the spin-Hamiltonian near the LAC, the 
electro-nuclear spin projection m = m s +mj is an approximately good quantum and the levels in figure [Sl7| 
have been labelled accordingly. In the region of the LAC, the transverse hyperfine terms A± mix levels of 



the same spin projection m, such that the levels with m = — | are mixed, but the levels with m — an( i 
— | are approximately unmixed. Strictly, the m = + ^ and — ? levels do undergo a small LAC (as indicated 



in figure S17) induced by the relatively small E terms. The E term and anisotropic hyperfine coupling also 
result in small admixtures of the | + 1, ±1/2) hyperfine levels to the eigenstates near the LAC. We proceed 
by formulating a basic model, where we neglect the E term and assume isotropic hyperfine coupling. In 
this case, the mixed states are composed from only two components | — 1, +1/2) and |0, —1/2) (total spin 
m = — |) and they can be quantified by the coefficients a and (3 (with \a\ 2 + |/3| 2 = 1) as denoted in figure 

Em 

We now use this hamiltonian to construct a basic 10-level rate-equation model that describes nuclear 
spin polarization by optical pumping and optical nuclear spin readout [51 We first neglect electronic 
spin relaxation and assume that the population rates of each of the triplet electronic spin states are equal. 



This remaining transitions are summarized in figure S18 To simplify the notation, we denote the triplet 
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Figure S18: 10-level rate equation model for optical nuclear spin polarization and readout. 



population rate far away from the LAC by 7 and the decay rates out of the electronic triplet states 70, 7-i, 
7+1. The population and decay rates of the unmixed states are 



\E,±l/2) -> 
|£,+l/2) 

|£?,-l/2)-> 

I + 1 ± 1/2} - 
I0.+1/2)- 

| - 1, -1/2} - 



+ l,±l/2) 
->|0,+l/2> 
-1,-1/2) 

> |G,±l/2) 

> \G,+l/2) 

> |G,-l/2> 
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7 

7+i 

7o 

7-i 



(S.17) 



The transitions rates for the mixed states \I) = a\-l, +l/2)+0|O, -1/2) and \II) = a|0, — 1/2}— /3| — 1, +1/2) 
are easily found by considering the matrix elements of the transition operator D. For example the decay 
rate from |J) to ground state +1/2 is found as 



7|/>|G,+i/2> = |<G, +l/2|D(a| - 1, +1/2) + /3|0, -1/2))| 2 

= \a(G, +l/2\D\ - 1, +1/2) + 0(G, +1/2\D\0, -1/2) | 2 
= \a{G, +l/2\D\ -1,+1/2)| 2 - |a|Vi 



(S.18) 
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Spin polarization proceeds in two steps. First, we perform optical pumping until steady state is reached. 
Thereafter the system is allowed to decay to the ground state. Before we evaluate the numerical solution to 
the rate-equations, we provide a qualitative argument and consider the case at the center of the LAC, where 
a = (3 = 2 -1 / 2 . Under this condition, the total population rates of all triplet hyperfine states are equal. 
Thus, the steady state population is determined by the decay rates. Specifically, population will be trapped 
predominantly in the longest lived hyperfine level. The lifetimes Tj = 1/7, of the triplet states are 

1 + 1, ±1/2): 7^ « 200ns 

a\ - 1, +1/2) + /?|0, -1/2) : (l/2 7o + l/2 1 _ l y 1 a UOOns 

|0,+l/2): 7,7 1 « 2500ns (S.20) 

a\0, -1/2) +1/2) : (l/2 7o + 1/2^-x)- 1 « 1400ns 

|-l,-l/2): jZl w 1000ns 

The | + l±l/2) states do not contribute imbalance of the decay rates for the +1/2 and —1/2 nuclear spin 
manifolds and can be ignored. The role of the mixed states is twofold. First, they provide nuclear spin flip 
transitions, and second, due to their larger decay rates to the +1/2 ground state (see above), they result 
in an overall increased rate from the —1/2 to the +1/2 manifold. The unmixed states provide nuclear spin 
polarization by their different lifetimes. Specifically, population will accumulate predominantly in the long 
lived state |0, +1/2). Both mechanisms act in the same direction, such that, in the steady state under strong 
optical pumping, the largest fraction of the population will be in |0, +1/2), a relatively large population will 
also be in the mixed states, minor population will be in the | — 1, +1/2) state and negligable population 
in the | + l,±l/2) states. During a subsequent decay (after the pump laser is turned off), all population 
from the |0, +1/2) will decay to the +1/2 ground state and the mixed states will preferentially decay into 
the +1/2 ground state. This will further increase the nuclear spin polarization. This process is described 



in terms of the change in the two ground state populations during a pulse sequence in figure S19 A full 
numerical analysis predicts a nuclear spin polarization (n +1 / 2 — n_i/2)/(n + i/2 + "--1/2) ~ 60%. 

We now consider nuclear spin readout. Starting from ground state +1/2 or ground state —1/2, we apply 
optical pumping until steady state is reached and evaluate the integrated fluorescence, that is given by the 
excited state populations (multiplied by the optical decay rate). The nuclear spin states will have different 
probability to be trapped in the long lived shelving state |0, +1/2) and thus have a different overall fluo- 
rescence intensity. The mixed states do not contribute to readout contrast, since their population rates are 
equal and as soon as they are populated, the nuclear spin information is lost. Likewise, the | + 1, ±1/2) 
states do not contribute to readout contrast and the only relevant states are the two unmixed states. Ap- 



preciating the lifetimes evaluated in Eq.(S.20), we expect that the nuclear spin \G, +1/2) will correspond to 
low fluorescence. A full numerical analysis predicts a fluorescence contrast between \G, +1/2) and \G, —1/2) 
of about 7%. The numerical analysis also predicts that the time to reach steady state is about 6 — 8/is, as 
observed in experiment. 

We now combine the results for nuclear spin polarization and readout contrast. The nuclear spin signal 
will be the product of both efficiencies. From this we expect a signal contrast of about 4%. This slightly 
overestimates the experimentally observed contrast of about 1.5 %. This likely reflects the approximations 
made in this basic model. Next, we can consider the sign of the signal contrast contrast. Experimentally, 
we observe a maximum of the nuclear Rabi oscillations at zero microwave pulse length, indicating that the 
nuclear spin state that we polarize in corresponds to bright fluorescence. By contrast, the basic model 
suggests the opposite case, where the polarized state (|0, +1/2)) corresponds to low fluorescence. This is 
likely also a result of the simplifications made in above basic model. 

Several aspects are likely to be relevant in a full description of the nuclear spin polarization- and readout 
contrast. One aspect is the small admixture from the | + 1 ± 1/2) states into the eigenstates at the LAC. Due 
to the comparably fast decay rate of the these states, the relatively small admixture might become relevant. 
A second aspect are possible small differences among the population rates. While these are not important 
for the electron spin polarization, here they could become relevant. After designing a precise measurement 
for the population rates far away from the LAC, these could be included in a complete rate model. Lastly, 
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Figure S19: Build-up of population in the ground state. Each curve represents the time-evolution of 
two nuclear spin state populations during a pulse cycle. Note that 60 % polarization is already achieved 
after the first cycle (not shown) and this plots shows a time-evolution after many cycles. 



electron-spin relaxation might be considered. Electron-spin relaxation would result in additional transition 
rates among the triplet hyperfinc levels. 
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